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Abstract 

We present the results of deep spectroscopy for the central region of the cluster lens SDSS J1004+4112 
with the Subaru telescope. A secure detection of an emission line of the faint blue stellar object (component 
E) near the center of the brightest cluster galaxy (Gl) confirms that it is the central fifth image of the 
lensed quasar system. In addition, we measure the stellar velocity dispersion of Gl to be cr* = 352 ± 13 
km s~^. We combine these results to obtain constraints on the mass Mbh of the putative black hole (BH) 
at the center of the inactive galaxy Gl, and hence on the Mbh-c* relation at the lens redshift zi — 0.68. 
From detailed mass modeling, we place an upper limit on the black hole mass, Mbh < 2.1 x IO^^Mq at 
la level (< 3.1 x IO^^Mq at 3a), which is consistent with black hole masses expected from the local and 
redshift-evolved Mbh-ct, relations, Mbh ^ 10^ - 10^°Mq. 
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1. Introduction 

Quasars lensed by foreground clusters of galaxies serve 
[ as a powerful probe of the mass distributions of clus- 
■ ters. SDSS J1004+4112 is the first example of such 
quasar-cluster lens systems (Inada et al. 2003; Oguri et 
al. 2004; Ota et al. 2006). It consists of four bright quasar 
images {z = 1.734) with a maximum image separation of 
14^6 produced by a massive cluster at z = 0.68. In ad- 
dition to the quasar images, multiply- imaged background 
galaxies are also observed (Sharon et al. 2005). An ad- 
vantage for using lensed quasars is that they can provide 
time delays, which have been actually detected for this 
system (Fohlmeister et al. 2007; Fohlmeister et al. 2008), 
containing unique information on the lens potential. 

What makes the quasar lens SDSS J1004-I-4112 partic- 
ularly unique is the probable central fifth image (com- 
ponent E) of the lensed quasar system. High-resolution 
Hubble Space Telescope (HST) Advanced Camera for 
Survey (ACS; Clampin 2000) images show a blue point 
source located 0."2 from the center of the brightest cluster 

* Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan, and 
observations (GO-9774, GO-10509, and GO-10716) made with 
the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. 



galaxy Gl (Inada et al. 2005). The central lensed image, 
if confirmed, has several important implications for the 
central structure of lensing objects. For instance, central 
images can constrain the inner mass profiles of lensing 
galaxies, particularly the masses of the central supermas- 
sive black holes (e.g., Mao et al. 2001; Rusin & Ma 2001; 
Keeton 2003; Rusin et al. 2005). While the correlation be- 
tween the black hole mass (Mbh) and the stellar velocity 
dispersions (a*) of the host galaxies has been established 
for local galaxies (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000; Tremaine et al. 2002), measurements of the 
redshift evolution for the correlation, derived for galaxies 
hosting active galactic nuclei, are still controversial (e.g., 
Peng et al. 2006; Shen et al. 2008; Woo et al. 2008). Thus 
independent constraints from central images are thought 
to provide insights into this relation. 

In this Letter, we present results of two deep spectro- 
scopic follow-up observations for SDSS J1004-I-4112, con- 
ducted at the Subaru 8.2-meter telescope. One is spec- 
troscopy of component E to confirm its lensing nature, 
and the other is spectroscopy of galaxy Gl to measure 
its velocity dispersion. The latter observation is par- 
ticularly important in separating the mass distribution 
of dark matter from that of baryons (e.g.. Sand et al. 
2008). We then attempt mass modeling of the lens sys- 
tem. Together with the measurement of a.^, for galaxy Gl, 
it provides a direct constraint on the Mbh-c* relation at 
the lens redshift of z = 0.68. In what follows, we assume 
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Fig. 1. HST ACS F435W image of SDSS J1004+4112 with 
the exposure time of 5400 sec, taken on 2005 December 
19 under the GO-10509 program (PI: C. S. Kochanek). 
Components labeled by A-D indicate the 4 lensed quasar im- 
ages, and "Gl core + E" indicates the central fifth image E 
and the center of galaxy Gl. The gray rectangle shows the slit 
direction for spectroscopy of the fifth image. The inset shows 
an expanded view of the central region of galaxy Gl. North 
is up and East is left, and the pixel scale is 0."05 pixel"^. See 
Inada et al. (2005) for HST images in other (redder) bands. 

a standard flat universe with = 0.26, fl\ = 0.74, and 
Ho = 72kms^^Mpc-i (e.g., Tegmark 2006). 

2. Observations and Data Analysis 

2.1. Fifth Image 

Wc conducted spectroscopic observations of the central 
fifth image (component E; see Figure 1) with the Faint 
Object Camera And Spectrograph (FOCAS; Kashikawa 
et al. 2002) at the Subaru 8.2-meter telescope on 2007 
January 22. We used the 300B grism, the L600 filter, 
and a I'.'O- width slit, under the 2x2 (spatial x spectral) 
on-chip binning mode. With this configuration the wave- 
length coverage is from 3650A to 6OOOA, which covers a 
strong emission line (C IV) of the lensed quasar with a 
spectral resolution of R ~ 400 and a pixel scale of 0"207 
pixcl^^ . Although component E is quite close to the cen- 
ter of galaxy Gl at 2; = 0.68 (0"2 from the center of Gl), 
our blue spectroscopy minimizes the contamination from 
Gl. The slit direction was aligned to pass through com- 
ponent D (~ lO'.'O from component E, see Figure 1) of 
SDSS J1004+4112 which we adopt as a reference for the 
quasar emission lines. Given the faintness {B ~ 24.5) of 
component E, we used a total exposure time of 16,200 sec 
taken in excellent seeing conditions (FWHM of ~ 0''6). 
After removing cosmic rays using the Lacos_spec task 
(van Dokkum 2001), we extracted 1-dimensional spectra 
of components E and D using the standard IRAF tasks. 
The spectra, normalized by the continuum level of each 
image, are shown in Figure 2. Wc find that the spectrum 
of component E shows a clear emission line just at the 
wavelength of the CiV emission line of component D. We 
note that there is no sky emission line near the wavelength 
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Fig. 2. Spectra of components E (black) and D (gray), nor- 
malized by each continuum level. In component E, an emis- 
sion line is detected at the same wavelength as the C IV emis- 
sion line of component D (also indicated by a vertical dotted 
line). The origin of the absorption features around 4200A and 
44OOA are unknown. 

of the C IV emission line (4200 A). Thus our spectroscopic 
observations unambiguously confirm that component E is 
the central fifth image of the lens system, representing 
the first spectroscopic confirmation of central odd images 
among lensed quasars. 

2.2. Velocity Dispersion of Gl 

To measure the velocity dispersion of galaxy Gl, we 
also used FOCAS but with a different setup. We adopted 
the 600-650 VPH grism (High 650), the SY47 filter, and a 
0"4-width slit. The wavelength coverage in this setup is 
from 53OOA to 77OOA, covering the Call H&K and G-band 
absorption lines of galaxy Gl at z = 0.68. The slit direc- 
tion was aligned along the major axis of galaxy Gl (160.8 
deg East of North). The spectrum, with a total exposure 
of 21,600 sec, was obtained on 2007 January 21, under 
the 3x1 on-chip binning mode (a spatial scale of 0''311 
pixel"^). The setup was chosen to achieve a spectral res- 
olution of R ~ 2500. The seeing was ~ 0."7 FWHM. Again, 
cosmic rays were removed by the Lacos_spec task and the 
1-dimensional spectrum of galaxy Gl was extracted by 
the standard IRAF tasks. 

We measure the velocity dispersion of Gl following the 
procedure described in Falco et al. (1997). We adopted 
HD 8491, HD 83085, HD 94247, and HD 126778 in the 
coude feed spectral library (Valdes et al. 2004) as template 
stars. We rebinned and smoothed the template spectra to 
match the observed spectral resolutions of our data. Both 
the spectrum of Gl and the resolution-matched template 
stars are normalized by the continuum levels. We then ap- 
ply the Fourier cross-correlation method (Tonry & Davis 
1979) of the IRAF FXCOR task, to the templates con- 
volved with Gaussians. The resultant calibration curves 
describe the relation between the input velocity disper- 
sion and the FWHM of the output cross-correlation peak 
measured by the FXCOR task. We used velocity disper- 
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sions from 200 kms^^ to 450 kms^^ with an interval of 
10 kms^^. The FXCOR task was also applied to the nor- 
malized rest-frame spectrum of Gl to obtain the width 
of the G-band absorption line^, which we convert to the 
velocity dispersion using linear interpolation of the cali- 
bration curves. We find that 4 template stars yield sim- 
ilar velocity dispersions for Gl of 338-367 kms^^. From 
the average and scatter (no weight for the 4 independent 
measurements), we determine the velocity dispersion to 
be tr* = 352 ± 13 kms^^. Our measurement is consistent 
with (T* expected from the observed Faber-Jackson rela- 
tion, ~ 320 ±40 kms-i (Liu et al. 2008). 

3. Constraints on the A^bh-c* Relation 

The location and the brightness of the confirmed cen- 
tral image allow us to place constraints on the mass of 
the supermassive BH hosted by Gl. This is made pos- 
sible by the fact that central point masses (de-)magnify 
or even suppress central images (e.g., Mao et al. 2001; 
Rusin et al. 2005). Indeed the technique was applied to 
three-image galaxy-scale lens PMN J1632— 0033 to derive 
an upper limit on A/bh in the lensing galaxy (Winn et al. 
2004), although the lack of the redshift measurement for 
the lensing galaxy makes the interpretation of this result 
somewhat difficult. 

Before conducting detailed mass modeling, we can esti- 
mate the BH mass range from the distance between com- 
ponent E and the center of Gl. It is given by the mass at 
which the Einstein radius of the central BH becomes com- 
parable to this distance of ~ 0"2, corresponding to a BH 
mass of Mbh 1.5 x IO^^Mq. Thus, we naively expect 
Mbh<10^°Mq. 

To obtain more quantitative constraints, we fit the 
quasar images with the following mass model. We adopt 
an elliptical version of the (Navarro et al. 1997; NEW) 
density profile for the dark matter distribution of the lens- 
ing cluster. Changing the inner slope of the dark matter 
component is expected to have little effect on the fifth 
image, because the center of the dark halo appears to be 
offset from Gl (Oguri et al. 2004). The central position, 
ellipticity, position angle, concentration parameter, and 
total mass of the dark matter distribution are treated as 
free parameters. The central galaxy Gl is assumed to be 
an isothermal ellipsoid with varying core radius. We con- 
strain the position angle to coincide with that of the light 
(161 ± 10 deg East of North), but with no prior on the 
ellipticity. In addition we added a Gaussian prior on the 
velocity dispersion of the isothermal ellipsoid Gl from our 
measurement, a = 352 ± 13 kms~^. The 16 galaxies other 
than Gl, which we identify as cluster members using their 
colors, are included as truncated isothermal ellipsoids with 
an exact scaling relation, cr* oc L^/"* and rtrunc oc L^^^. The 
overall normalization of rtrunc is left as a free parameter. 
We take the ellipticities and position angles for these mem- 
ber galaxies from measurements in the HST ACS images. 

^ The G-band is the only line which we can use, because the Ca II 
H&K absorption lines are not suitable for measuring velocity 
dispersions (Tonry 1998; Ohyama et al. 2002). 
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Fig. 3. The critical curves of our best-fit mass model. The 
observed positions of five quasar images are marked by filled 
squares, which are reproduced by our model nearly perfectly. 
The best-fit parameters are MjMpw = 2 X 10^^ Mq, cj^fw = 2, 
a^NFW = 6"184, j/NFW = 6"639, CNp-yv = 0.3, 9nfw = —11°, 
eoi = 0.71, 6»Gi = 152°, ctqi = SSSkms"!, rcoro.Gi = 0"l, 
7 = 0.21, = 82°, and rtrunc(L = ici) = 3."2. 

Finally, we allow for an external shear following Oguri et 
al. (2004) to achieve better fits. For observational con- 
straints, we adopt the positions and position errors of five 
quasar images and Gl derived in Inada et al. (2005). We 
also include flux ratios of quasar images as constraints, but 
with large errors of cr(mx — wa) = 0.3 (X represents im- 
ages B-D) and a{mE — fnp^) = 0.8, in order to allow small 
flux ratio changes due to time delays and microlensing. 
In addition, observed time delays, AtAB = 40.6 ± 1.8 days 
(Fohlmeister et al. 2007) and AiAC = 821.6 ±2.1 days 
(Fohlmeister et al. 2008), are included. 

We find that our model reproduces the data quite well 
(see Figure 3 for the critical curves of our best-fit model). 
With no BH, the best-fit chi-square is = 2.57 for 3 de- 
gree of freedom (20 constraints and 17 parameters). Thus 
the central BH is not required to fit the fifth image. We 
then place the BH (modeled by a point mass) at the cen- 
ter of galaxy Gl, and we re-perform minimizations for 
each Mbh- We note that the central BH can produce 
a sixth image near the BH. In what follows we ignore 
the sixth image as it is much fainter than the observed 
fifth image in most of the situations considered here. We 
derive the upper limit of A/bh from the differences, 
= X^(-^bh) — Xmin- Figurc 4 plots Ax^ as a function 
of Mbk , from which we derive constraints on the BH mass 
of A/bh < 2.1 x 10^° AT© at la and Mbh < 3.1 x 10^°Mq 
at 3cr. The steep increase of for large A/bh means that 
the mass of the BH is limited to being a small fraction of 
the stellar mass interior to image E. 

As shown in Figure 4, our constraints on A/bh are con- 
sistent with the expected local A/bh-c* relation (Tremaine 
et al. 2002) and also those inferred from the extrapolation 
of a possible redshift evolution measured by Woo et al. 
(2008). Although the constraints are model-dependent, 
in the sense that they are derived assuming a specific 
parametric mass model, we note that our mass model 
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Fig. 4. The difference (Ax^) as a function of the BH mass 
A^BH- The steep rise of Ax^ provides upper limits on Mbh- 
The shaded regions indicate the values of Mbh expected from 
the measured stellar velocity dispersion cr* = 352 ± 13 km s~^: 
The "local" value is estimated from the Mbh-o"» relation of 
Tremaine ct al. (2002), A/bh = 10**-^^(o-,/200kms-l)* 02A4'Q 
with an intrinsic dispersion of ~ 0.25 dex, whereas the "evo- 
lution" value is computed in the same way as "local" , ex- 
cept for an additional offset of the BH mass due to a 
possible redshift evolution measured by Woo et al. (2008), 
AlogMBH =3.11og(l + 2) = 0.70 for z = 0.68. 

is quite flexible near the center of Gl because the core 
radius is a free parameter. Indeed, the best-fit core ra- 
dius is correlated with Mbh; the core radius evolves from 
- 0."1 for A/bh = to ~ 0."3 for Mbh = 3 x 10^°Mq. The 
rough agreement of this constraint with a simple order- 
of-magnitude argument also suggests that an upper limit 
to the BH mass of approximately 2 x IO^^Mq is generic. 
However, improved constraints may be obtained by adding 
constraints on the sixth image which we ignored, because 
the sixth image becomes brighter with increasing A/bh- 
For instance, the flux of the sixth image is ~ 20% of E in 
our best-fit models for A/bh ~ 2 x 10^° Af©, whereas the 
HST image suggests that it cannot be ^ 10% of E. 

4. Summary 

We have presented the results of two spectroscopic ob- 
servations at the Subaru telescope. With the first observa- 
tion, wc confirmed that the central point source found by 
Inada et al. (2005) is indeed the central fifth image of the 
lensed quasar by detecting the C IV emission line of the 
quasar. This represents the first spectroscopic confirma- 
tion of the central odd image. In the second spectroscopic 
observation, we determined the stellar velocity dispersion 
of galaxy Gl to be cr* — 352 ± 13 kms~^. 

We used these results to derive constraints on the A/bh- 
(7* relation at the lens redshift z ~ 0.68, assuming the 
presence of a BH at the center of the brightest cluster 
galaxy Gl. With a parametric model which successfully 
reproduces the model constraints, we obtained limits of 
A/bh < 2.1 x IO^^Mq at la and Mbh < 3.1 x IO^^Mq at 
3(7 that are consistent with A/bh expected from the ex- 
trapolation of the known A/bh-c* relation. It is worth 



noting that these constraints arc derived for an "inactive" 
galaxy with no nuclear activity. Current studies of the 
redshift evolution of the A/bh-c* relation make use of the 
Balmer line widths, and therefore are restricted to BHs in 
active galaxies. Given possible differences in the A/bh-c* 
between active and inactive local galaxies (Greene & Ho 
2006), constraints on A/bh for distant inactive galaxies 
from gravitational lensing will be essential to fully under- 
stand the origin of the relation. 
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